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_ . I. INTRODUCTION 
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pL| ■ The adiabatic theorem QLQIEQ is one of the oldest and most widely used general tools in quantum mechanics. 
' The theorem concerns the evolution of systems subject to slowly varying Hamiltonians. Roughly, its content is that if 
a state is an instantaneous eigenstate of a sufficiently slowly varying H at one time then it will remain an eigenstate at 
later times, while its eigenenergy evolves continuously. When the slowness assumption is relaxed transitions become 
' weakly allowed 0, S S 13 • The role of the adiabatic theorem in the study of slowly varying quantum mechanical 
systems spans a vast array of fields an d ap plications, such as the Landau-Zener theory of energy level crossings in 
molecules 0,0], quantum field theory ^IJ, and Berry's phase [T^. In recent years geometric phases 01 have been 
proposed to perform quantum information processin g |14| |l5l Il6j| . with adiabaticity assumed in a number of schemes 
for geometric quantum computation (e.g., [TtI ITsl llfll l2flj). Additional interest in adiabatic processes has arisen 
in connection with the concept of adiabatic quantum computing, in which slowly varying Hamiltonians appear as 
a promising mechanism for the design of new quantum alg orithms and even as an alternative to the conventional 
quantum circuit model of quantum computation [2l|, [2^, l23| . 

More recently, in Ref. |24j . the adiabatic theorem was generalized to the case of open quantum systems, i.e., 
quantum systems coupled to an external environment. Instead of making use of eigenstates of the Hamiltonian, 
adiabaticity is defined through the Jordan canonical form of the generator of the master equation governing the 
dynamics of the system. This new framework allowed for the derivation of an adiabatic approximation which includes 
the case of systems evolving in the presence of noise. This issue is particularly important in the context of quantum 
information processing, where environment induced decoherence is viewed as a fundamental obstacle on the path to 
^ \ the construction of quantum computers (e.g., |25j]L 

The aim of this paper is to review the adiabatic approximation in quantum mechanics for both closed and open 
quantum systems as well as to point out how an incorrect manipulation of the adiabatic theorem can yield an 
inconsistent result. Indeed, in a recent paper entitled "Inconsistency in the application of the adiabatic theorem" j2(J 
the authors argue that there may be an inconsistency in the adiabatic theorem for closed quantum systems. We show 
here how this inconsistency can be resolved. Related discussions can be found in Refs. [271 l28l |29| . 



II. THE QUANTUM ADIABATIC APPROXIMATION FOR CLOSED SYSTEMS 



A. Condition on the Hamiltonian 



Let us begin by reviewing the adiabatic approximation in closed quantum systems, which evolve unitarily through 
a time-dependent Schrodinger equation 

H(t)\m)=i\m), a) 

where H(t) denotes the Hamiltonian and \ip{t)) is a quantum state in a D-dimensional Hilbert space. We use units 
where Ti — 1. For simplicity we assume that the spectrum of H(t) is entirely discrete and nondegenerate. Thus we 
can define an instantaneous basis of eigenenergies by 



H(t)\n{t)) =E n (t)\n(t)), 



(2) 
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with the set of eigenvectors \n(t)) chosen to be orthonormal. In this simplest case, where to each energy level there 
corresponds a unique eigenstate we can define adiabaticity as the regime associated to an independent evolution of the 
instantaneous eigenvectors of H(t). 

This means that instantaneous eigenstates at one time evolve continuously to the corresponding eigenstates at later 
times, and that their corresponding eigenenergies do not cross. In particular, if the system begins its evolution in a 
particular eigenstate |n(0)) then it will evolve to the instantaneous eigenstate \n(t)} at a later time t, without any 
transition to other energy levels. 

It is conceptually useful to point out that the relationship between slowly varying Hamiltonians and adiabatic 
behavior can be demonstrated directly from a simple manipulation of the Schrodinger equation: recall that H (t) can 
be diagonalized by a unitary similarity transformation 

H d (t) = U- 1 (t)H(t)U(t), (3) 

where Hd(t) denotes the diagonalized Hamiltonian and U(t) is a unitary transformation. Multiplying Eq. by 
U~ l {t) and using Eq. J2J we obtain 

H^d^i^d-iU- 1 ^), (4) 

where \ip)d = U^ 1 ^) is the state of the system in the basis of eigenvectors of H(t). Upon considering that H(t) 
changes slowly in time, i.e. dH(t)/dt w 0, we may also assume that the unitary transformation U(t) and its inverse 
C/ _1 (t) are slowly varying operators, yielding 

Hd(t)m))d = i\m) d . (5) 

Thus, since Hd(t) is diagonal, the system evolves separately in each energy sector, ensuring the validity of the adiabatic 
approximation. 
Now let 



g nk (t) = E n (t) - E k (t) 



(0) 



be the energy gap between level n and k and let T be the total evolution time. One may then state a general validity 
condition for adiabatic behavior as follows: 



max 

0<t<T 



(k\H\n) 



9nk 



< min \gnk\- 

0<t<T 



(7) 



Note that the left-hand side of Eq. Q has dimensions of frequency and hence must compared to the relevant physical 
frequency scale, which can be proved to be given by the gap g n k 0,H3|. [In fact, Eq. J7J will be seen to be a direct 
consequence of the adiabatic condition derived in Subsection lllBI ] The interpretation of the adiabaticity condition {7J) 
is that for all pairs of energy levels, the expectation value of the time-rate-of-change of the Hamiltonian, in units of 
the gap, must be small compared to the gap. For a discussion of the adiabatic regime when there is no gap in the 
energy spectrum see Ref. |3l| . 

In order to obtain Eq. (0, let us expand \ip(t)) in terms of the basis of instantaneous eigenvectors of H (t): 



(8) 



n = l 



with a n (t) being complex functions of time. Substitution of Eq. JSJ into Eq. Q and multiplying the result by (k(t)\, 
we have 



a-k 



-J2^n(k\h)e- l Io dt ' 9 ^ t ' ) . 



(9) 



A useful expression for {k\h), for k ^ n, can be found by taking a time derivative of Eq. (j2J) and multiplying the 
resulting expression by (k\, which reads 



<fc|n) = min) (»/fc). 



9nk 



(10) 
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Therefore Eq. can be written as 

(k\H\n) -i f* dt' g nh (t') 
a k = -a k (k\k) - 2_^dn - e Jo H " fcW . (11) 

Adiabatic evolution is ensured if the coefficients a k {t) evolve independently from each other, i.e., if their dynamical 
equations do not couple. As is apparent from Eq. (|ll(l . this requirement is fulfilled when the condition (JTJ is imposed. 

In the case of a degenerate spectrum of H(t), Eq. 1)11) [) holds only for eigenstates \k) and \n) for which E n ^ E k . 
Taking into account this modification in Eq. {TJ| , it is not difficult to see that the adiabatic approximation generalizes 
to the statement that each degenerate eigenspace of H(t), instead of individual eigenvectors, has independent evolution, 
whose validity conditions given by Eq. JJJ are to be considered over eigenvectors with distinct energies. Thus, in general 
one can define adiabatic dynamics of closed quantum systems as follows: 

Definition II. 1 A closed quantum system is said to undergo adiabatic dynamics if its Hilbert space can be decomposed 
into decoupled Schrodinger-eigenspaces with distinct, time- continuous, and non-crossing instantaneous eigenvalues of 
H(t). 



B. Condition on the total evolution time 



A very useful alternative is to express the adiabaticity condition in terms of the total evolution time T. We shall 
consider for simplicity a nondegenerate H(t); the generalization to the degenerate case is also possible. Taking the 
initial state as the eigenvector |m(0)), with a m (0) = 1, the condition for adiabatic evolution can be stated as follows: 



T > Jj, (12) 



where 



T = n 1 ?^ K fc ( s )l ~7~~ l m ( g ))l ' & = \9mk{s)\ ■ (13) 

o<s<i as o<s<i 

Eq. I|12l) can be interpreted as stating the total evolution time must be much larger than the norm of the time- 
derivative of the Hamiltonian divided by the square of the energy gap. It gives an important validity condition for the 
adiabatic approximation, which has been used, e.g., to determine the running time required by adiabatic quantum 
algorithms pHl22ll23]| . By using the time variable transformation l|ltj[l . one can show that Eq. (|12fl is indeed equivalent 
to the adiabatic condition J7J) on the Hamiltonian. 

To derive Eq. i fEfy , let us rewrite Eq. l(TT|) as follows : 



4 («*(*) e-*»M) (14) 
where Jk{t) denotes the Berry's phase 0] associated to the state \k): 

7fc(*)=i / dt'(k(t')\k(t')). (15) 



Now let us define a normalized time s through the variable transformation 

t = sT, < s < 1. (16) 
Then, by performing the change t — > s in Eq. 114|) and integrating we obtain 

o fc 00e-*»W = Ofc(0 ) - V f^s'M^e-^i:'^^'')^ (17) 

^Jo 9nk{s') 



where 



F nk (s) = a n (s) (fc( a )|^W|n(s)> e^W. (18) 
ds 
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However, for an adiabatic evolution as defined above, the coefficients a n (s) evolve without any mixing, which means 
that a n (s) fa a n (0) e^"W. Therefore 

F nk (s) = a n (0) (k(s)\^^-\n(s)) e -^«-T««). (19) 
as 

In order to arrive at a condition on T it is useful to separate out the fast oscillatory part from Eq. (|17|l . Thus, the 
integrand in Eq. I|17|l can be rewritten as 

Fnk(s') -iT J° ds"g nk (s") _ i 
9nk{s') T 

Substitution of Eq. lj2U|) into Eq. (|T7|) results in 

a k (s)e-^ = a k (0) + ^ ~ F ° ^ + f ds ' ^ "'^'^¥Ss) ' ^ 

n^k \ g nk\ V ) 9 n k\ S ) JO " S 9 n k\ S ) / 

A condition for the adiabatic regime can be obtained from Eq. <|21|l if the last integral vanishes for large T. Let us 
assume that, as T — > oo, the energy difference remains nonvanishing. We further assume that d{F nk (s')/ g^k(s')} /ds' 
is integrable on the interval [0, s]. Then it follows from the Riemann-Lebesgue lemma [33j that the last integral in 
Eq. H21fl vanishes in the limit T —> oo (due to the fast oscillation of the integrand) [34j . What is left are therefore 
only the first two terms in the sum over n ^ k of Eq. I|21l) . Thus, a general estimate of the time rate at which the 
adiabatic regime is approached can be expressed by 

9 

where 

F= max \a n (0) (k(s)\^W-\n(s)}\, g= min \g nk (s)\ , (23) 

0<s<l ds 0<s<l 

with max and min taken over all k and n. Eq. (|12|l is then obtained as the special case when the system starts its 
evolution in a particular eigenstate of H(t). 



d ( F n k(s') -iT J 3 ds"g nk (s") 

ds' M C ° 



-tT f 3 ds"g nk (s")d_ ( F nk {s') 

ds'\g 2 nl .(s') 



(20) 



C. Higher-order corrections to the adiabatic approximation 



When the Hamiltonian of a quantum system changes slowly, but not extremely slowly, the degenerate eigenspaces 
of H(t) (or individual eigenvectors in the case of nondegenerate spectrum) will not evolve completely independently 
from each other and, therefore, the dynamical equation Q will weakly couple distinct eigenspaces of H(t). Then, for 
non-extremely slowly varying Hamiltonians, the adiabatic solution is actually a zeroth-order approximation and higher- 
order corrections must be considered. Some higher-order adiabatic approximation methods have been proposed 0, 
IE • Here we shall review, for the non-degenerate case, the method proposed by Wu in Ref. Q . Let us begin by 
expanding the state vector \ip(t)} in the instantaneous eigenbasis, as in Eq. JSJ. Then, by using the normalized time 
s introduced in Eq. (|16|l we obtain the following matrix form for the Schrodingcr equation 



diP(s) 
ds 

where K(s) is an anti-Hermitian matrix with elements 



K(sW(s), (24) 



K mn (s) = - (m(s)| |n(s))exp (iT J ds' g mn (s'^j . (25) 

The matrix K(s) can be separated into a diagonal matrix D(s) and an off-diagonal matrix O(s), yielding 

K{s) = D{s) + 0(s). (26) 
The evolution operator U (s) for the system satisfies the equation 

= K(s)U(s), (with U(0) = 1) (27) 

ds 
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which, after integration and use of Eq. (|26|l . becomes 

U(s) = l+ [ ds 1 [D(s 1 )+0(s 1 )]+ f d Sl f 1 ds 2 [D( Sl ) + 0( Sl )}[D{s 2 ) + 0{s 2 )}+ ... (28) 
Jo Jo Jo 

Now let us define 

C/ (0) (s) = 1+ f d Sl D( Sl )+ [ dsx [ ds 2 D( Sl )D(s 2 ) + (29) 
Jo Jo Jo 

which involves only the diagonal parts. Moreover, for n > 0, we denote U^ n \s) as the sum of all the integrals with 
n off-diagonal O(s) factors in the integrand. Therefore the evolution operator can be expanded in powers of the 
off-diagonal matrices O(s) as U(s) — J2 n =o U^^s). It can be shown || that the n th term U^ n \s) can be expressed 
through the lower order term U^ n ~^(s) by means of the recurrence equation 

U {n) {s)= f ds'C/(°)(s')0(s')C/ (n " 1) (s')- (30) 
Jo 

The expression above means that, by knowing the zeroth-order evolution operator U^°'(s), which exactly yields the 
adiabatic approximation, one can obtain [/W(s), and then U^- 2 \s) and so on. The adiabatic case corresponds to 
no transitions, while a correction £A n )(s) of order n > 1 implies the existence of n transitions between different 
energy levels |8|- From this perspective one can interpret the adiabatic approximation as the zeroth order term in a 
perturbation theory in the number of transitions between energy levels connected by the time-varying Hamiltonian. 



III. THE QUANTUM ADIABATIC APPROXIMATION FOR OPEN QUANTUM SYSTEMS 

In this section we review our recently introduced generalization of the adiabatic theorem to the case of open 
quantum systems |2 ll . The motivations for considering such a generalization are many. The most fundamental is that 
the concept of a closed system is, of course, an idealization, and in reality all experimentally accessible systems are 
open. Thus applications of the adiabatic theorem for open systems include, among others, geometric phases (where 
open system effects have received considerable recent attention, e.g., Refs. H3), quantum information processing, 
and molecular dynamics in condensed phases. 

In the following we first introduce notation for open systems, then discuss the generalized adiabatic theorem. 



A. The dynamics of open quantum system 

Consider a quantum system S coupled to an environment, or bath B (with respective Hilbert spaces 7is,7is), 
evolving unitarily under the total system-bath Hamiltonian H$ b ■ The exact system dynamics is given by tracing over 
the bath degrees of freedom [3(j 

p(t) = Tr B [U(t) PsB (0)UHt)}, (31) 

where p(t) is the system state, psb(0) = p(0) ® ps(0) is the initially uncorrelated system-bath state, and U(t) = 
Texp(— i H,gR(t')dt') (T denotes time-ordering). Such an evolution is completely positive and trace preserving 
|36l l37l l38| . Under certain approximations, it is possible to convert Eq. (13 1|) into the convolutionless form 

p{t) = C{t)p{t). (32) 

An important example is 

P (t) = -i[if(t),p(t)] + i^([r i (t),^)rt(t)] + [r,(t) P (i) ) rl( i )]). (33) 

i=l 

Here H(t) is the time-dependent effective Hamiltonian of the open system and Ti(t) are time-dependent operators 
describing the system-bath interaction. In the literature, Eq. (|33() with time-independent operators I\ is usually 
referred to as the Markovian dynamical semigroup, or Lindblad equation [3^ . I38ll39l l4fl | [see also Ref. for a simple 
derivation of Eq. (1331) from Eq. I|31|l ]. However, the case with time-dependent coefficients is also permissible under cer- 
tain restrictions |42| . The Lindblad equation requires the assumption of a Markovian bath with vanishing correlation 



6 



time. Equation (|32() can be more general; for example, it applies to the case of non-Markovian convolutionless master 
equations studied in Ref. Here we will consider the class of convolutionless master equations <|32l) . In a slight 

abuse of nomenclature, we will henceforth refer to the time-dependent generator C(t) as the Lindblad superoperator, 
and the as Lindblad operators. 

Conceptually, the difficulty in the transition of an adiabatic approximation from closed to open quantum systems 
is that the notion of Hamiltonian eigenstates is lost, since the Lindblad superoperator - the generalization of the 
Hamiltonian - cannot in general be diagonalized. It is then not a priori clear what should take the place of the 
adiabatic eigenstates. This difficulty was solved in Ref. 24] by introducing the idea that this role is played by 
adiabatic Jordan blocks of the Lindblad superoperator. The Jordan canonical form ^4| , with its associated left- and 
right-eigenvectors, is in this context the natural generalization of the diagonalization of the Hamiltonian. In this 
direction, it is convenient to work in the superoperator formalism, wherein the density matrix is represented by a 
D 2 -dimensional "coherence vector" 

\p)) = (pi Pi ••• Pz>» )*> (34) 

and the Lindblad superoperator C becomes a D 2 x D 2 -dimensional supermatrix [2^. We use the double bracket 
notation to indicate that we are not working in the standard Hilbert space of state vectors. More generally, coherence 
vectors live in Hilbert-Schmidt space: a state space of linear operators endowed with an inner product that can be 
defined, for general vectors u and v, as 

(u,v) = ((u\v))=^Tr(u^v). (35) 

where TV is a normalization factor. Adjoint elements ((v | in the dual state space are given by row vectors defined 
as the transpose conjugate of \v)): ((v\ = {v\,V2, ■■■,v* D - 2 ). A density matrix can then be expressed as a discrete 
superposition of states over a complete basis in this vector space, with appropriate constraints on the coefficients so 
that the requirements of Hermiticity, positive semi-definiteness and unit trace of p are observed. Thus, representing 
the density operator in general as a coherence vector, we can rewrite Eq. 132(1 in a superoperator language as 

C(t)\p(t))) = \p(t))), (36) 

where C is now a supermatrix. This master equation generates non-unitary evolution, since C(t) is non-Hermitian 
and hence generally non-diagonalizable. However, one can always transform C into the Jordan canonical form |44j |. 
where it has a block-diagonal structure. This is achieved via the similarity transformation 

Cj(t) =S- 1 (t)C(t)S(t), (37) 

where £j{t) = diag(Ji, J m ) denotes the Jordan form of C(t). The Jordan blocks J Q , of dimension n a , are always 
of the form: 



Jot 



f\ a 1 ... \ 
A Q 1 ... 



... A a 1 

V A Q J 



(38) 



To each Jordan block are associated a left and a right eigenvector with eigenvalue A Q , which can in general be 
complex. The number m of Jordan blocks is given by the number of linearly independent eigenstates of £(i), with 
each eigenstate associated to a different block J a . Since C(t) is non-Hermitian, we generally do not have a basis 
of eigenstates, whence some care is required in order to find a basis for describing the density operator. It can be 

shown |24[ that instantaneous right ||P^(i)))| and left j((£a^(i)| j bases in the state space of linear operators can 

always be systematically constructed, with the following suitable features: 

• Orthonormality condition: 

((£W(t)\vV\t))) = j{{£®\Df))j = 8 afl 6V. (39) 

• Invariance of the Jordan blocks under the action of the Lindblad super-operator: 

C{t) |2#>(t)» = \T>U-»(t))) + A a (t) |2#5(t)». (40) 
«£W(t)| C(t) = + «^(*)| A a (t). (41) 

with \T>a = and ({Ea 1 "^ = 0. Here the subscripts enumerate Jordan blocks (a G {1, ...,m}), while the 
superscripts enumerate basis states inside a given Jordan block G {0, n a — 1}). 
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B. Adiabatic conditions for open quantum system 



Before stating explicitly the conditions for adiabatic evolution, we provide a formal definition of adiabaticity for 
the case of open systems: 

Definition III.l An open quantum system is said to undergo adiabatic dynamics if its Hilbert- Schmidt space can 
be decomposed into decoupled Lindblad-Jordan-eigenspaces with distinct, time- continuous, and non-crossing instanta- 
neous eigenvalues of C(t). 

This definition is a natural extension for open systems of the idea of adiabatic behavior. Indeed, in this case 
the master equation l|32f> can be decomposed into sectors with different and separately evolving Lindblad-Jordan 
eigenvalues. The more familiar notion of closed-system adiabaticity is obtained as a special case when the Lindblad 
superoperator is Hermitian: in that case it can be diagonalized and the Jordan blocks all become one-dimensional, 
with corresponding real eigenvalues (that correspond to energy differences). The splitting into Jordan blocks of the 
Lindblad superoperator is achieved through the choice of a basis which preserves the Jordan block structure as, for 

example, the sets of right j|Z^'(i))} j and left j((£i^(t)| j vectors introduced above. Such a basis generalizes the 

notion of Schrodinger-eigenvectors. Based on this concept of adiabaticity, we state below (without the proofs) several 
theorems which have been derived in Ref. Y2M. 



Theorem III. 2 A sufficient condition for open quantum system adiabatic dynamics as given in DeRnition UII. 1\ is: 

« 1, (42) 



max 

0<s<l 




(& +P ~ 1] \&t Sp) )) 



where s = t/T is the scaled time and 



w^ a (t) = A„(t) - X a (t), 5, = Vfc s (with S = 0), III V I = >>••>;> ( 43 ) 




and where Xp \ a , with i and j denoting arbitrary indices associated to the Jordan blocks a and 8, respectively. 

The role of the energy differences that appear in the equations for the closed case is played here by the (in general 
complex- valued) difference between Jordan eigenvalues top a , while the norm of the time-derivative of the Hamiltonian 
is replaced here by the norm of 4^, evaluated over and inside Jordan blocks. Eq. (|42|l gives a means to estimate the 
accuracy of the adiabatic approximation via the computation of the time derivative of the Lindblad superoperator 
acting on right and left vectors. The norm used in Eq. 1|42[) can be simplified by considering only the term with 
maximum absolute value, which results in: 



Corollary III. 3 A sufficient condition for open quantum system adiabatic dynamics is 

« 1, (44) 



■n a np 

y 0<s<l 



hf':T ,J max 



«^ +p - x) lf \vt Sv) )) 



P+S p 



where the max is taken for any a ^ 0, and over all possible values of i £ {0, n a — 1}, j 6 {0, np — 1}, and p, 
with 

(l+j)!(n a -i)! [0> 

The factor A^°"^ given in Eq. I|45|l is just the number of terms of the sums in Eq. I|42|) . We have included a 
superscript np, even though there is no explicit dependence on np, since j £ {0, ...,np — 1}. 

Furthermore, an adiabatic condition for a slowly varying Lindblad super-operator can directly be obtained from 
Eq. lO, yielding: 

Corollary III. 4 A simple sufficient condition for open quantum system adiabatic dynamics is C « 0. 



8 



Note that this condition is in a sense too strong, since it need not be the case that L is small in general (i.e., for all 
its matrix elements). 

Just as in the closed-systems case, one can also express the condition for adiabaticity in terms of the total time of 
evolution. To this end, we expand the density matrix for an arbitrary time t in the instantaneous right eigenbasis 



as 



m nff—X 



W))) 



1 p rt 



(46) 



P=i j=o 



where m is the number of Jordan blocks and n@ is the dimension of the block J p. We emphasize that we are assuming 
that there are no eigenvalue crossings in the spectrum of the Lindblad superoperator during the evolution. We also 
define 



v^ p \s)= P f{s){{e^{s)\ 



djC(s) 
ds 



\V ( t Sp \s))) 



(47) 



and 



JW*) = / U/3<x{t')dt'. 



(48) 



Then the following adiabatic time condition can be established [24J: 

Theorem III. 5 Consider an open quantum system governed by a Lindblad superoperator C(s). Then the adiabatic 
dynamics in the interval < s < 1 occurs if and only if the following time conditions, obtained for each coefficient 
Pa\s), are satisfied: 



T 3> max 

0<s<l 



/3|A^#A Q j,p 



p+S p +l 



V^ p) (s)e TQ ^ 



(0) 



P+Sp+l 
J l3a 



(s) 



+ / ds'e TQ ^ 



d V^(s') 



ds' ,,,p+ s p+ 1 



J f3a 



(S>) 



(49) 



Theorem IIII. 51 provides a very general condition for adiabaticity in open quantum systems. Equation (|49|l simplifies 
in a number of situations. 

• Adiabaticity is guaranteed whenever (s) vanishes for all X a =/= Xp . 

• Adiabaticity is similarly guaranteed whenever Vg^ (s), which can depend on T through pi , vanishes for all 
A Q , Xp such that Re(Qp a ) > and does not grow faster, as a function of T, than exp(T| Rettp a \) for all A a , A^ 
such that Re(0 / a Q( ) < 0. 

• When Re(fl/3 a ) = and Im(rip ) ^ the integral in inequality Q49JI vanishes in the infinite time limit due to 
the Riemann-Lebesgue lemma [33j], as in the closed case discussed before. In this case, again, adiabaticity is 
guaranteed provided Pp\s) [and hence Vg^(s)] does not diverge as a function of T in the limit T — > oo. 

• When Re(fi / g Q ) > 0, the adiabatic regime can still be reached for large T provided that Pp (s) contains a 
decaying exponential which compensates for the growing exponential due to Re(O i g a ). 

• Even if there is an overall growing exponential in inequality Ij49(l , adiabaticity could take place over a finite time 
interval [0, T*] and, afterwards, disappear. In this case, which would be an exclusive feature of open systems, 
the crossover time T» would be determined by an inequality of the type T ^> a + 6exp(cT), with c > 0. The 
coefficients a, b and c are functions of the system-bath interaction. Whether the latter inequality can be solved 
clearly depends on the values of a, 6, c, so that a conclusion about adiabaticity in this case is model dependent. 

A simpler sufficient condition can be derived from Eq. I)49|l by considering the term with maximum absolute value 
in the sum. This procedure leads to the following corollary: 

Corollary III. 6 A sufficient time condition for the adiabatic regime of an open quantum system governed by a 
Lindblad superoperator C(t) is 



T » 



max 

0<s<l 



< 3P) (Q) 

p+Sp + l 



" /3a 



(0) 



p+S p + l 
J 0o 



is) 



, y(ijp)( A 



ds' , .p+s p +i 

pa 



(s') 



(50) 
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where max is taken over all possible values of the indices X a Xp, i, j, and p, with 

mZ-"= E E E II E M-^Pk^w (51 » 

where Ap a denotes the number of Jordan blocks such that X a ^ A/3. 

Further discussion of the physical significance of these adiabaticity conditions, as well as an illustrative example, 
can be found in Ref. 24] . The application of the adiabatic theorem for open quantum systems to problems in 
quantum information processing (e.g., in the context of adiabatic quantum computing |2lL l22l |23Tp and geometric 
phases |45l 1461 l47j . seems particularly appealing as a venue for future research. 



IV. THE MARZLIN-S ANDERS INCONSISTENCY 



The adiabatic theorem can be deceptively simple when it is not carefully interpreted. In a recent paper entitled 
K.-P. Marzlin and B.C. Sanders argue that there may be an inconsistency in the adiabatic theorem for closed quantum 
systems |26l . when the change in instantaneous adiabatic eigenstates is significant. Here we simplify their argument 
and show where exactly is the fallacy that leads one to conclude that there is such an inconsistency. 



A. The condition for adiabaticity revisited 

Let us consider a quantum system evolving unitarily under the Schrodinger equation Q . At the initial time to the 
system is assumed to be in the particular instantaneous energy eigenstate \Eo(to)}. For a general time t the evolution 
of the system is described by 

\iP(t)) = U(t,t )\E (t )}, (52) 

where U(t,t a ) — Texp(— i f t H{t')dt') is the unitary evolution operator. Assuming that the Hamiltonian changes 
slowly in time and that \Eo(to)) is non-degenerate, the adiabatic theorem implies that 

\m)=e- i ^ Eo e i ^\E (t)) (53) 



where / Eq = f tg Eo(t')dt' and the Berry's phase (3q is given by /3q = i J{Eq\Eo). Therefore the substitution of 
Eq. into the instantaneous eigenbasis of H(t), defined by H{t)\E n {t)) — E n {t)\E n (t)) , yields 

H(t)m)) = E (t)\iP(t)}, (54) 

which simply states that the wave function is an instantaneous eigenstate of the Hamiltonian in the adiabatic regime. 
Substituting Eq. I|54|l into the Schrodinger equation one obtains 

i\i]j)=E (t)\m)- (55) 

It is important to observe that the above equation has been derived by using the fact that the adiabatic solution must 
be an instantaneous eigenstate of the Hamiltonian. Moreover one can see that the adiabatic wave function is really a 
solution of the adiabatic Schrodinger equation by substituting Eq. (|53|l into Eq. (|55|l . from which one obtains 

(1 - |S (*)X^o(*)I) |^o(*)> = 0. (56) 

In order to show that Eq. I|56|l is obeyed in the adiabatic regime we can project this equation by multiplying it by 
each instantaneous basis vector {E n it) \ : 

(E n (t)\ (1 - \E (t)}(E (t)\) \E (t)) = 0. (57) 

Therefore we obtain that the above equation is satisfied, for each n, if the adiabatic constraints [see Eq. Q] are 
obeyed 



(E n (t)\E (t)) 



E (t)-E n (t) 



< 1, (n ^ 0) (58) 
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B. The inconsistent step 

Now suppose that we wish to solve the adiabatic Schrodinger equation l|55(l , but (incorrectly) ignore the fact that it 
has been derived by assuming that \"ip(t)} is an instantaneous eigenstate of H{t). Then, imposing the initial condition 
l^(*o)) — one easily finds that Eq. Q55f) is satisfied by the following wave function: 

m))=e- i f tEo \E (t )). (59) 

However, this \ip(t)) is clearly inconsistent with Eq. Q54JI and therefore is an illegal solution, since it generally is not 
an instantaneous eigenstate of the Hamiltonian. Indeed, the general adiabatic solution is given by Eq. (|53[) . which 
includes the Berry's phase and \E (t)), as opposed to |i^o (^o)) - In fact, if we take Eq. ijSlH) as the adiabatic wave 
function and substitute it into Eq. I|53|) we obtain 

e i ^\E (t)) = \E (to)}. (60) 

Then multiplying Eq. (JHUI by (E (t )\ 

e i0 °(E Q (t Q )\E o (t)) = l. (61) 

This inconsistency is precisely the one claimed by Marzlin and Sanders in Ref. [26^. Eq. (6). Note that this result 
is obtained there in a somewhat more complicated manner, by considering the adiabatic solution of a time-reversed 
wave function |V>(£)) = U'(t, t )\E n (t )). We note that the solution for \rjj(t)) [their Eq. (4)] is very similar to our 
Eq. (|59|l . Hence, in the same way that we have been led to an inconsistent result due to a deliberately wrong adiabatic 
solution for |V>(*))i Ref. [2(| has been led to an inconsistent solution for their |V>(i)}. 

V. CONCLUSION 

We have reviewed the adiabatic dynamics of both closed and open quantum systems. In the case of closed systems the 
adiabatic limit is the case where initial Schrodinger-eigenspaces evolve independently, without any transitions between 
eigenspaces; this limit can be relaxed and a perturbation theory can be developed in the number of transitions. In the 
case of open systems the notion of Schrodinger-eigenspaces is replaced by independently evolving Lindblad-Jordan 
blocks. A corresponding perturbation th eory has not yet been developed. We have also shown that the inconsistency 
in the adiabatic theorem claimed in Ref. 26] is a consequence of an improper adiabatic solution for the wave function. 
One arrives at an inconsistent result by taking the instantaneous adiabatic eigenstates and integrating them over all 
time using the adiabatic Schrodinger equation. The adiabatic theorem remains a valuable and consistent tool for 
studying the dynamics of slowly evolving quantum systems. 
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